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ABSTRACT 

We have made continuum and spectral line observations of several outflow 
sources in the Mon OBI dark cloud (NGC 2264) using the Heinrich Hertz Tele- 
scope (HHT) and ARC 12m millimeter-wave telescope. This study explores the 
kinematics and outflow energetics of the young stellar systems observed and as- 
sesses the impact star formation is having on the surrounding cloud environment. 
Our data set incorporates ^^CO (3-2), ^^CO (3-2), and ^^CO (1-0) observations 
of outflows associated with the sources IRAS 06382+1017 and IRAS 06381+1039, 
known as IRAS 25 and 27, respectively, in the northern cloud complex. Comple- 
mentary 870 fim continuum maps were made with the HHT 19 channel bolometer 
array. Our results indicate there is a weak < 0.5 % coupling between outflow ki- 
netic energy and turbulent energy of the cloud. An analysis of the energy balance 
in the IRAS 25 and 27 cores suggests they are maintaining their dynamical in- 
tegrity except where outflowing material directly interacts with the core, such as 
along the outflow axes. 

Subject headings: stars: formation, ISM: individual - Mon OBI (NGC 2264) - 
Clouds - ISM: jets and outflows - ISM: kinematics and dynamics - ISM: molecules 
- radio lines: ISM 



1. Introduction 

The Mon OBI dark cloud lies toward the galactic anticenter and has been the subject 
of several large-scale molecular hne surveys (Blitz 1979), (Oliver, Masheder, & Thaddeus 
1996). Due to its location, there is little source confusion from intervening foreground and 
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background clouds. The relative proximity of this complex, 800 pc according to the work of 
Walker (1956), has also contributed to its survey appeal. The region of Mon OBI associated 
with NGC 2264 has itself been the subject of several studies, including an unbiased ^^CO (1 — 
0) survey and the discovery of nine molecular outflows (Margulis & Lada 1985), (Margulis, 
Lada, & Snell 1988). Margulis, Lada, & Young (1989) presented the results of an IRAS 
survey of 30 discrete far-infrared sources, including IRAS 25 and IRAS 27. The IRAS Point 
Source Catalog (PSC) identifies IRAS 25 and IRAS 27 with objects IRAS 06382+1017 and 
IRAS 06381+1039, respectively. A systematic search for dense gas with CS observations 
revealed that these objects were associated with a ~ 500 — 700 clump (Wolf-Chase, 
Walker, & Lada 1995). 

Wolf-Chase et al. (2003) used HIRES processing of IRAS data, ^^CO (2-1) observations, 
and SCUBA continuum images to identify individual cores embedded within NGC 2264 and 
concluded that IRAS 25 and IRAS 27, initially thought to be Class I protostellar objects 
(Margulis, Lada, & Young 1989), consist of multiple sources in different evolutionary stages. 
They also find evidence for multiple outfiow driving sources in the ^^C0(2 — 1) data. In 
a wide field imaging survey of the NGC 2264 region of Mon OBI in the infrared, H^, and 
"^^CO (3 — 2), Reipurth et al. (2004) report several new Herbig-Haro (HH) objects near the 
star forming cores found by Wolf-Chase et al. (2003) and attempt to pinpoint their origins. 

A primary goal of this work is to determine physical characteristics and outfiow prop- 
erties of IRAS 25 and 27. With this information, we hope to gain a better understanding of 
the star formation history, dynamics, and impact that these protostellar objects have on the 
surrounding cloud environment. In Section (2.), we describe the observations presented in 
this work. Section (3.) addresses the results of the continuum and spectral line observations 
of each of the sources, IRAS 25 and 27. In Section (4.), we discuss velocity centroid studies of 
the molecular line data and draw energetics comparisons to assess the impact star formation 
is having upon the cloud complex. A summary of our findings is presented in (5.) 



2. Observations 

2.1. Continuum mapping with the HHT 

The continuum data presented in this work were obtained with the facility 19 channel 
array operating at the HHT on Mt. Graham, Arizona. This array, developed by E. Kreysa 
and collaborators, contains 19 individual broadband bolometric detectors sensitive around a 
central wavelength of A = 870 fira (345 GHz) and arranged in two concentric hexagonal rings 
surrounding a central pixel. The field of view of the array is ~ 200" in diameter. Observations 
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were made during 2003 January 12 - 14*'^ and May 6 - 8*'*, and 2004 April 21 - 22"'^ using 
On The Fly (OTF) mapping. With this technique, fully sampled 450" x 300" maps were 
made surrounding the sources IRAS 25 and 27 in NGC 2264. A total of six maps of IRAS 27 
and three maps of IRAS 25 were combined to produce the data sets presented here. During 
observations, the pointing was checked each hour by observing a planet, such as Venus, 
Saturn, and Jupiter. For all of the data, the pointing varied by less than 4". These maps were 
corrected for the effects of atmospheric opacity using "skydip" scans taken prior to beginning 
and upon completion of each map. This correction is accomplished by fitting a polynomial 
function to the optical depth, at 870 //m as a function of elevation. The bolometer data 
were flux calibrated with images of Mars, Saturn, and the known submillimeter secondary 
calibrator source IRAS 16293-2422 (Sandcll 1994) taken using a continuum on/off sequence 
technique. We estimate the flux calibration to be accurate within the ~ 20% level. 

2.2. Spectral line observations 
2.2.1. HHT data 

In February and May of 2004, we obtained ^^CO (3 — 2) and ^^CO (3 — 2) spectral line 
data with the HHT toward known outflow regions, NGC 2264 O and NGC 2264 H, associated 
with IRAS 25 and 27, respectively. A summary of important observational parameters for 
these data is presented in Table 1. A set of 300" x 300" ^^CO (3 - 2) fully sampled OTF 
maps were made surrounding the two IRAS sources. The pointing was monitored similar 
to the continuum observations and varied < 5". A total of two separate maps, each with 
two polarization data, were made for NGC 2264 O &: H. After careful calibration, maps 
were combined to produce the data described in Table 1. In addition, ^^CO (3 — 2) Absolute 
Position Switched (APS) spectra were taken toward the locations of the IRAS sources as 
published in Margulis, Lada, & Young (1989). A total integration time of 140 s was spent on 
source for each APS ^^CO observation. The 2a RMS of all observations has been included 
in Table 1. 

The HHT facility dual-polarization, double-sideband 345 GHz receiver was used with two 
Acousto-Optical Spectrometer backends, each possessing 1 GHz total bandwidth and 2048 
channels. Individual APS spectra were taken toward the positions of IRAS 25 and 27 at the 
beginning of each observing session and used to cross cahbrate data sets taken on different 
days. Orion A (0:1950 = 05'* 32"^ 47.0^ ^1950 = -05° 24' 21") was used to calibrate the dual 
polarization data before they were combined. For ease of comparing data sets from different 
telescopes, the main beam efficiency {rj^b) of the 345 GHz receiver was measured and used 
to convert antenna temperature corrected for atmospheric extinction, T4, to a main beam 



-4- 



brightness temperature, T^b, according to T^b — TX/rimb- Measurements of Saturn were 
used to find rjmb using Equation 1, where T^lplanet) is the planet's brightness temperature, 
D is the angular diameter of the planet when it was observed, and 6 is the telescope beam 
size. Measurements of Saturn made during observations resulted in rjmb values that varied 
~20%. Based on this variation, we adopt a value of rjmb ~ 60% for the 345 GHz receiver 
during this time. 

TU planet) ,^ /--D^, ^m~^ 

Vmb = IT 1 - exp(—^ln2)] 1 

Tji{planet) 0^ 

2.2.2. 12m data 

APS ^^CO (1—0) observations of IRAS 25 and 27 were made at the ARO 12m millimeter- 
wave telescope near Tucson, AZ during April 2002. A summary of the observational infor- 
mation is presented in Table 1. A total integration time of 120 s was spent on source for 
each observation and 2a RMS values are listed in Table 1. The 12m dual-polarization, 
single-sideband 90 — 116 GHz receiver and Millimeter Autocorrelator (MAC) backend with 
300 kHz total useable bandwidth and 8192 channels were used. Similar to the HHT data, 
both polarizations were combined in order to improve the signal-to-noise using Orion A as a 
cross-calibrator. The 12m makes measurements in terms of T^, the observed antenna tem- 
perature corrected for atmospheric attenuation, radiative loss, and rearward scattering and 
spillover. In order to make meaningful comparisons between data sets, these measurements 
were converted to a T^b scale using a procedure documented in Appendix C of the 12m 
manual. According to the manual, Tmb — T^/1.08 for the 90 — 116 GHz receiver at 115 GHz. 
With this conversion, the 12m data were compared with HHT spectra after the data sets 
were convolved to the same beam size. Measurements of T^b for the 12m data are accurate 
to within ~ 20%. 



3. Results 

3.1. Continuum Analysis 

Results of the 870 /xm HHT continuum data analysis are presented for IRAS 25 & 27. 
Table 2 lists the fluxes found at the positions of sources IRAS 25 & 27 and the associated 
continuum cores. The cores, labeled with SI, S2, and S3 extensions, were initially idcntifled 
as emission peaks in 850 /^.m and 450 /xm SCUBA maps (Wolf-Chase et al. 2003). We derived 
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870 /im peak fluxes within multiple apertures for each of the sources, including 45" (for easy 
comparison with Wolf-Chase et al. (2003)) and 3'x5', the size of the IRAS error eUipse at 
100 urn. 

Combining our 870 //m continuum data with bolometric observations of other observers 
and IRAS, we construct FIR SEDs (see Figure 1) for dust cores with 870 //m fluxes presented 
in Table 2. When IRAS data were used, the fluxes were taken from the PSC At 60 /im & 
100 fim, the IRAS the beam sizes are larger than the IRAS 25 source, indicating fluxes of the 
entire clump emission were found. When constructing the IRAS 25 FIR SED, measurements 
of the 870 fim emission of the whole clump were used. In order to incorporate SCUBA 
450 /xm & 850 /xm data, we estimated the total flux found at these wavelengths by summing 
pubhshed fluxes within 45" apertures of the 25 SI and S2 cores. For IRAS 27, 60 jj,m. & 
100 iim IRAS fluxes were incorporated in the FIR SED with the 870 iim fluxes measured 
within 45" apertures around the continuum cores IRAS 27 (27 SI), S2, and S3. Published 
SCUBA 450 fim & 850 /im measurements made within the same 45" apertures were included. 
When available (for sources 27 SI & S2) 1300 fim fluxes were incorporated. A summary of 
continuum measurements for each source at all wavelengths where data are available is in 
Table 3 along with associated references. 

In order to probe physical conditions of the continuum cores, FIR SEDs were flt with 
expressions of the following form (Hildebrand 1983): 

= nsB,{T){l - e^-') (2) 

where is the source size and 

BAT) - ( j^tyfi) & r. = (./.„)'' (3) 

In Equation 3, Biy{T) is the Planck function, r,^ describes the frequency dependent optical 
depth, P is the spectral index governing dust grain emissivity properties, and Uo is the 
frequency where Tjy = 1. 

The results of fltting flux values to Equations 2 and 3 are illustrated in Figure 1 and 
Table 4 summarizes the details numerically. We have only attempted to flt the cold dust 
component (A > 60 /xm) of the continuum emission that is responsible for the majority of 
the flux at mm and submm wavelengths. In all cases, except for IRAS 25, each continuum 
core was well flt by a single temperature blackbody (BB), similar to Equation 3. In order 
to obtain a good least-squares flt to each of these dust cores (IRAS 27 (27 SI), 27 S2, and 
27 S3), source size (Qs), temperature (T), and spectral index (/3) parameters were varied 
iteratively until a best flt solution was found. 
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Fitting SED curves to IRAS 25 flux measurements proved more complex, and a two 
temperature component BB function was used. In order to obtain a good, physically plausible 
fit to these data, both dust temperature components, T, and j3 indices were selected and 
held constant while the source sizes, Qg, were iteratively fit. The goodness of fit did not 
strongly depend upon the values of (3 chosen for the two components and therefore the 
resulting fit does not reflect unique or precise values of spectral index. Since spectral index 
is probably constant within the same molecular cloud complex, our choice of (3 used to fit 
the IRAS 25 emission was guided by solutions from iterative fits of IRAS 27 (27 Si), 27 S2, 
and 27 S3 SEDs. The temperatures of the two components fit to IRAS 25, 10 K and 40 K, 
were selected to be similar to the ~ 30 K obtained from using Wien's displacement law to 
estimate the characteristic dust temperature for IRAS 25. The two component BB fit was 
found to be rather sensitive to the values of temperature chosen. Based upon this sensitivity, 
the temperatures of the components are probably accurate to within 5K. 

Luminosities (Lfir) of continuum cores were found using Lpjji = 3.1 x 10~^'^D'^ f Fudu 
(Walker, Adams, & Lada 1990) where Lfir is in L©, D is in pc, F^, is in Jy, and u is in GHz. 
Dust filling factors, fd = Qg/^m, were estimated by comparing values of Qg derived from 
the SEDs to half-power 870 fim source sizes, Qm- Table 4 presents the results of the FIR 
SED model fits and parameters for each of the cores examined. Using techniques similar 
to Hildebrand (1983), physical properties of the dust cores, including dust optical depth at 
870 //m {ti, — T87o), column density (Nh^), and total gas mass {MH2) were derived. Number 
density values (n^jj) were found using the relation Mjj2/(4/37r i^^mijj) where is the 
mass of a hydrogen molecule and R is the radius over which M^^ was found. The visual 
extinction was estimated using the relation = 1.1 x 1Q~'^^Nh2 (Bohlin, Savage, & Drake 
1978). The results of this analysis are presented in Tabic 5. In making these computations, 
we assumed a gas to dust ratio of 100 in mass, a grain density of p = 3 g cm^"^, and a 
dust particle radius of a = 0.1 /im. These values arc consistent with those determined by 
Hildebrand and arc commonly used for molecular clouds in the galaxy. 

At 870 fim, Figure 1 shows that the majority of the fiux observed for IRAS 25 is due to 
the colder 10 K BB component. It is reasonable to obtain an upper limit to the amount of 
mass in the IRAS 25 core region by assuming that all fiux at these wavelengths is attributable 
to the cold gas. Using the analysis tools described above, a limit of 140 Mq was obtained for 
the IRAS 25 core (Table 5). Although this mass limit is large compared with the IRAS 27 
core, it is not unexpected. The observed flux at 870 //m for IRAS 27 (Figure 1) is well 
described by a 30 K BB source and is weaker by about a factor of two when compared with 
IRAS 25 870 //m fluxes (Table 3). At these wavelengths, we observe comparatively more flux 
from a colder source, implying larger column density and core mass values for IRAS 25. 
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3.2. Spectral Line Observations 

In this section, the results of -"^^CO (3 — 2) & ^"^CO (3 — 2) analyses are presented for 
IRAS 25 & 27. The C0(3 — 2) spectra for both objects are single peaked features and possess 
broad line wings (~40 km s~^ in extent for IRAS 25 and ~50 km s^^ for IRAS 27). In 
order to investigate the spatial distribution of the outflow molecular gas, emission line wings 
were defined by comparing the ^^CO (3 — 2) & ^^CO (3 — 2) spectra (see Figure 2) and then 
used to make contour maps. The ^^CO (3 — 2) line wings set the maximum extent of low 
velocity emission associated with the outflows. High velocity wings extended to where the 
^^CO (3 — 2) emission merged with the noise. Figure 2 illustrates these regions, including 
the line core, low, and high velocity wings for both IRAS 25 and 27 and Table 6 summarizes 
this information. In this table, the spectral line emission of the entire (full) red and blue 
wings (F RW & F BW) has been divided into low and high velocity (LV & HV) regions. 

3.2.1. IRAS 25 (IRAS 06382+1017) 

Contour maps of F RW and F BW integrated intensity were constructed for IRAS 25 
using If rw = £ fcm7-i ^^^^ k Ip bw = I^ioTZis-^'^A ^"^dv, respectively. Figure 3 
shows red (thin, solid contours) and blue wing (thin, dashed lines) ^^CO (3 — 2) emission 
associated with outflow NGC 2264 O surrounding IRAS 25. Contours are superimposed 
on an 870 /xm continuum OTF map (see Section (3.1)) and the figure includes locations of 
other sources discovered within the same region. Figure 3 shows two prominent blue outfiow 
lobes, SEE & WB, and an elongated red emission lobe. Although a single red outfiow 
lobe is detected, observations indicate that the emission is unresolved and two separate red 
lobes, ER & NWR, exist. Figure 4 is a contour map of the high velocity (HV) red and 
blue ^^CO (3 — 2) emission associated with NGC 2264 O. A large fraction of the eastern red 
emission (ER) disappeared for the high velocity gas, leaving a clear view of the compact 
northwestern red (NWR) emission. This indicates that the ER and NWR emission possess 
different velocity components and suggests that the emission belongs to two different outflow 
lobes. 

Our data support the existence of two ^^CO (3 — 2) outflows in the region. Figure 4 
aided in determining the orientations of the two outflows associated with NGC 2264 O. The 
high velocity gas traces an outflow oriented southeast (SEB) to northwest (NWR). The WB 
lobe and a portion of ER emission do not appear in this figure. This indicates these lobes 
have similar velocities and suggests they describe a second outfiow. The existence of multiple 
molecular outfiows within this region has been suggested by other authors (Piche, Howard, 
& Pipher 1995), (Wolf-Chase et al. 2003). Our observations of two ^^CO (3 - 2) outfiows fit 
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well with the work of Reipurth et al. (2004) and lend support to observations of recurring 
outflow activity taking place in the region surrounding IRAS 25 (Ogura (1995), Wolf-Chase 
et al. (2003), Reipurth et al. (2004)). 

For clarity, the information contained in Figure 3 is summerized in the model sketch 
of the IRAS 25 region shown in Figure 5. The location of IRAS 25 is denoted with a box, 
and continuum peaks 25 SI & 25 S2 associated with 850 ;um&450//m SCUBA emission are 
shown with "X" symbols. IRAS 25 has been identified as the driving source for a large scale 
(0.48 pc) Hcrbig-Haro fiow, HH 124 (Walsh, Ogura, & Reipurth 1992). Six knots of emission 
were observed (asterisks) and lie along the WB— ER outfiow axis. This region also includes 
an IR reflection nebula (thick, gray contour in Figure 5) associated with IRAS 06382+1017 
(Piche, Howard, & Pipher 1995) and two VLA sources (triangles) within the IRAS error 
ellipse (Rodriguez & Reipurth 1998). 

In this complex environment, it is difficult to determine sources of the ^^CO (3 — 2) 
outflows. Our observations and the work of Reipurth et al. (2004) indicate that the WB— ER 
outflow axis is well aligned with HH 124, suggesting that the IRAS object may be the source 
of this outflow emission. The axis of the SEB— NWR outflow is in proximity to many objects, 
including IRAS 25 and the submillimeter continuum sources. The source of this molecular 
outflow has not been determined (Reipurth et al. 2004), although it is oriented along the 
opening angle of the IR reflection nebula associated with IRAS 25 (see Figure 5). It is 
possible that two ^^CO (3 — 2) molecular outfiows emanate from IRAS 25; the WB— ER fiow 
oriented along HH 124, and a SEB— NWR outfiow along the IR refiection nebula opening 
angle. 

The SEB— NWR outfiow axis does not intersect well with IRAS 25. If the IRAS object 
is the source of the fiow, this may imply that the outfiow is redirected as it encounters dense 
ambient material. Figure 3 shows that the IRAS source is located near the edge of a dense, 
extended 870 /xm continuum core. We estimate how dense this material would have to be in 
order to redirect the outflow and make comparisons with the results of the continuum analysis 
{nH2 in Table 5). To redirect the flow, the clump pressure, given by the ideal gas law {rici k T) 
in LTE, must exceed the pressure of the outflowing material, Pof/{'^ R'^)- For an order of 
magnitude estimate, it is assumed the outflow acts mostly over an area encompassing its 1/2 
power contours. In this analysis, k is Boltzmann's constant, T = 10 K is the core temperature 
derived from continuum observations (see Table 4), Pqf is the outflow momentum rate (see 
Section (3.2.3)), and R is the average radius of 1/2 power outflow contours. The minimum 
core density capable of redirecting the outflow is Ud ~ lO^cm"^. This is large compared 
with (see Table 5), suggesting the density of the core is most likely not sufficient to 
influence outflows in the region. 
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It was suggested that IRAS 25 is a barely resolved binary system with a possible compo- 
nent ~ 1400 AU (1.8" at the assumed distance of NGC 2264) from the primary at a position 
angle of ~ 155° (Piche, Howard, & Pipher 1995). The presence of two outflows is consistant 
with IRAS 25 being a binary system and may explain the unusual morphology. The 870 fim 
continuum emission (Section (3.1)) can be used to estimate the mass of stellar objects in the 
system. Assuming that the bolometric luminosity {Lfir) of this object (see Table 4) equals 
the accretion luminosity {Lace — GMMacc/ R), an upper limit to the mass of the central 
object is obtained. We assumed a typical radius of it! = 2 x 10^^ cm, ~ 3 R©, for the hy- 
drostatic core of an accreting protostellar object (Stabler, Shu, & Taam 1980) and the mass 
infall rate of an isothermal cloud. Mace — 0.975 c^/G (Shu 1977), where Cg — yJWTJm^^ is 
the cloud sound speed. Since material is accreted from a reservoir in close proximity to the 
protostar, the warmer continuum component T = 40K (see Table 4) was used to calculate 
Cg- This simple analysis estimates 0.5 for the stellar object(s) associated with IRAS 25. 

3. 2. 2. IRAS 27 (IRAS 06381 +1039) 

Similar to (Section (3.2.1)), ^^CO (3 — 2) line wing integrated intensities were found for 
NGC 2264 H outflows surrounding IRAS 27 using Ip rw = /qT fe^T-i"' ^^dv k Ip bw ^ 
I-2oo"kLs-^'^A^'^dv. Figure 3 shows the outflow emission (F RW = solid, thin lines and 
F BW = thin, dashed contours) superimposed on a map of 870 ;um continuum emission. 
Other sources in the region are shown, including the location of the IRAS source (box) and 
870 & 450 //m SCUBA emission peaks, 27 SI, 27 S2, & 27 S3. Our observations indicate 
there is evidence for two ^^CO (3 — 2) outflows located near IRAS 27, supporting the recent 
work of Reipurth et al. (2004). Figure 4 shows the high velocity outflow emission that was 
used to help determine outflow orientations. The flgure indicates HV gas traces an outflow 
oriented southeast (SER) to northwest (NWB). The northeastern (NEB) and southwestern 
(SWR) lobes vanish in Figure 4 indicating that they have similar velocities and describe a 
second outflow. 

Figure 5 is a model sketch of the IRAS 27 system, showing the orientations of the 
^^CO (3 — 2) outflows. H2, Ha, and [SII] images of the region reveal an IR source surrounded 
by a reflection nebula (thick, gray contour) and two giant HH outflows, HH 576 and 577 
emanating from a core containing IRAS 27 (Reipurth et al. 2004). Figures 3 and 5 show 
submillimeter cores located along a dense ridge of 870 fim emission with molecular outflows 
oriented at large angles (nearly perpendicular) to the continuum ridge. It is difficult to 
determine the source of the ^^CO (3 — 2) outflows, given their proximity to the continuum 
cores and IRAS object. It is possible that both outflows originate from IRAS 27. 
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The axis of the NEB-SWR flow is not as well centered on IRAS 27 as the SER-NWB 
outflow. An order of magnitude calculation, similar to the analysis in Section (3.2.1), reveals 
the volume density of the environment, n/^j, is less than the density required to redirect the 
outflow {rici ~ 10^ cm~^). The alignment of HH and ^^CO (3 — 2) outflows prompted sug- 
gestions that IRAS 27 is a binary source with each component driving an outflow (Reipurth 
et al. 2004). If so, then following the same procedure as for IRAS 25, we estimate ~ 1 M© 
for the mass of the stellar object (s) associated with IRAS 27. 



3.2.3. Outflow Analysis 

Optical Depth, Excitation Temperature, & H2 Column Density 

Outflow characteristics and energetics analyses were carried out under assumptions of con- 
stant excitation throughout the mapped regions and Local Thermodynamic Equilibrium 
(LTE). With pointed observations toward IRAS 25 & 27 in ^^CO (3-2), ^^CO (3-2), and 
^^CO (1 — 0), measurements of optical depth and excitation temperature were made and H2 
column density was derived. Optical depth, r, was found using ^^CO (3 — 2) and ^^CO (3 — 2) 
observations. Assuming ^^CO & ^^CO have comparable excitation temperatures due to their 
similarity in rotational energy level structures, the ratio of their optical depths is the abun- 
dance ratio, Ta (see Equation 4). 

h2co 1 - e-^»^32 



13CO l-e-'^32^" 

Optical depth of the ^^CO was determined by making a ratio of the integrated intensities 
of the ^^CO & ^^CO spectra, I^co and hscOi respectively. An abundance ratio Va — 60 for 
^^CO to ^^CO was used (Kulesa 2003), Equation 4 was solved numerically. 

After convolving ^^CO (3 — 2) and ^^CO (1 — 0) data sets to the same beamsize (~ 1') 
and converting each temperature scale into Tmb as described in Section (2.2), line ratios of 
the spectra were used to determine T^^ (Walker 1991). In this analysis, the ^^CO (3 — 2) and 
^^CO (1—0) emission was assumed to originate from the same volume and since our data sets 
were convolved to the same beamsize, the emission fllling factors were the same. With this 
method, values of 20 K and 30 K were found for T^.^ toward IRAS 25 and 27. respectively. 

With estimates of optical depth and excitation temperature, the ^'^CO (3 — 2) column 
density was determined with the following expression (Groppi et al. 2004; Walker et al. 1988): 

6A: 2Z + 1 e^hv/{2kT,,) 
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In the above equation, I is the lower rotational state, u describes the upper state, k and 
h are respectively Boltzmann's and Planck's constants, u is the ^^CO (3 — 2) frequency, 
Av is the velocity interval over which the column density was calculated, and n is the 
molecular dipole moment expressed in esu. The above expression describes how to calculate 
the column density toward a position where the values of r & T^x are known. Pointed 
^^CO (3 — 2) observations were made toward individual IRAS sources, though ^^CO (3 — 2) 
data were not available over the regions mapped in ^^CO (3 — 2). For this reason, excitation 
temperature, optical depth, and column density values were obtained at the positions of the 
IRAS sources, although calculating this information directly over the mapped regions was 
not possible. Table 7 presents LTE properties of IRAS 25 and 27 calculated for red and blue 
line wings. 

In order to obtain column density information and ultimately mass and energetics infor- 
mation for outflows NGC 2264 O & H an emission-weighted area technique was used (Lada 
1985). The analysis assumed the observed intensity ratio of -"^^CO (3 — 2) and ^^CO (3 — 2) 
emission at the positions of the IRAS sources was constant over the regions mapped in 
^^CO (3 — 2). It was also assumed that T^x obtained at the positions of IRAS sources was 
constant over the mapped areas. The intensities of the ^^CO (3 — 2) emission measured 
toward the IRAS sources, where N^,thin was calculated, were compared with ^^CO (3 — 2) 
emission over the entire maps and used to weight the contribution of individual map positions 
to the total emission weighted area of the outflow, A^^ight- This procedure is described in 
Equation 6, where Aheam is the area of the telescope beam at 345 GHz, Ij.^y is the ^^CO (3 — 2) 
intensity at any map position, and Iiras 25/27 is the ^^CO (3 — 2) intensity at IRAS 25 or 
27, depending on the map examined. 



map 

^weight 



map 

-E-^^ean.( , ) (6) 



I- IRAS 25/27 
x,y I 



Mqf — — ^^-^^ Ayjeight fthin (7) 

rarH2 

Using the emission-weighted area method, the gas mass associated with the outflow was 
obtained using Equation 7. In this equation, ruHz is the molecular hydrogen mass and fthin is 
the ^'^CO (3 — 2) gas filling factor. The filling factor was calculated using Tex, the integrated 
intensity and optical depth of the ^^CO (3 — 2), and assuming LTE conditions. Filling factors 
obtained toward IRAS sources appear in Table 7. Since ^^CO (3 — 2) measurements were 
made toward single positions, fthin was assumed to be constant over the maps. In Equation 7, 
Va is the ^^CO/^^CO abundance ratio. We adopted = 60 in agreement with the work of 
Kulesa (2003). The CO/H2 abundance, r^jj = 1 x 10~^, is representative of the ISM and 
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many nearby molecular clouds. Results of our mass analysis for the red and blue NGC 2264 O 
& H outflow lobes are provided in Table 8. 

Outflow Energetics 

With outflow masses, we were able to derive other physical properties of NGC 2264 O 
& H, including momentum {Pqf = Mqf Vmax), kinetic energy {Eqf = 1/2 Mqf 'i^max)^ 
mechanical luminosity {Lqf = EoF/tdyn), mass outflow rate {Mqf = MoF/tdyn), and force 
{PoF = PoF/tdyn)- In this analysis, tdyn = Ling/vmax is the outflow dynamical age, is 
the outflow lobe physical extent measured from emission half-power contours, and v^ax is 
the characteristic velocity of the outflowing material. LTE outflow energetics depend largely 
upon the choice of Vmax which is subject to uncertainty due to projection effects of the 
outflow with respect to the line of sight. Determining a characteristic velocity representative 
of the bulk of the outflowing gas is therefore difficult. For this reason, lower and upper limits 
of the outflow energetics were determined using outer line wings of ^^CO (3 — 2) (LV RW 
& LV BW) and ^^CO (3 - 2) (F RW & F BW) for Vmax- The results of these calculations 
appear in Table 8. 

These energetics results are largely consistent with those of Margulis, Lada, & Snell 
(1988). Upper limit kinetic energy, momenta, and mechanical luminosities of the outflows 
are comparable to those typically associated with young high mass protostars, possessing 
Lfjr = 10^ - 10^ Lq and Mqf = tens of (Zhang et al. 2005). Since neither IRAS 25 
nor 27 is believed to be a high mass star forming object, these results are consistant with 
the interpretation that the outflows arise from multiple objects. 

4. Analysis &; Discussion 
4.1. Velocity Centroid Analysis 

In order to explore the effect the NGC 2264 O & H molecular outflows have on the cloud 

environment, we constructed velocity centroid maps of the ^^CO (3 — 2) gas. The velocity 
centroid is the value that equally divides the integrated intensity of a line profile. For this 
analysis, the velocity centroid was computed over the line core excluding emission in the 
wings (sec Table 6). A map of velocity centroid as a function of position was generated 
for NGC 2264 O & H (see Figure 6 a and b). The centroid velocity values have been 
represented with labeled, thin contours. Half power contours of the red and blue wing 
^^CO (3 — 2) emission are depicted with thick solid and dashed fines, respectively. There is 
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a small gradient in velocity centroid values in each map, Ai>cent = (7.65 — 6.6) km s ^ 1.1 
km s"^ and Avcent = (8.35 - 7.4) km s"^ fa 1 km s"^ for NGC 2264 O & H, respectively. 
There is no noticeable sign of rotation perpendicular to outflow axes to suggest the signature 
of an accretion disk. The centroid velocity gradients extend from the northeast corner of the 
plots to the southwest, roughly following the orientations of the outflow lobes. This suggests 
that NGC 2264 O & H outflows may dominate the energetics of their surroundings. 

In order to illustrate this quantitatively, the force of bulk motions of cloud material was 
calculated and compared with the force of the molecular outflow, Pqf- For this analysis, 
the mass of the cloud was computed over the line core for each ^^CO (3 — 2) map, Mc, using 
the emission weighted area approach described in Section (3.2.3). This mass calculation 
assumes the ^^CO (3 — 2) emission is optically thin over the line core. In order to assess the 
impact this might have on the analysis, the ^^CO (3 — 2) optical depth was examined as a 
function of velocity toward the position of the IRAS sources, where direct calculations of 
optical depth and column density were made. For IRAS 25. the -"^^CO (3 — 2) emission is 
marginally thick in the line core, with optical depth ~ 1 — 1.1. The mass associated with line 
core velocities, and therefore the corresponding force calculations, should serve as a lower 
limits for NGC 2264 O. Over the line core region, the ^^CO (3 — 2) emission toward IRAS 27 
is thin, with optical depth values ~ 0.3 — 0.35, so the mass associated with the line core, 
and therefore the corresponding force calculations, are most likely valid. 

Calculations of the force of bulk cloud motions were made using = Avcent/Tdyn, 
where Mc is the mass in the line core region, Avcent is the observed centroid velocity gradient 
along the outflow axis, and Tdyn is the dynamical age of the outflow. Using dynamical 
ages of the red and blue outflow lobes resulted in Pc = 6.6 — 11 x 10~^ yr""^ km 
for NGC 2264 O. This is more than ~5 times smaller than the force of the outflowing 
material, Pqf = 60 x lO^^Mgyr"^ km s~^, obtained by summing the contributions of the 
red and blue wings shown in Table 8. These comparisons were also made for energetics and 
dynamical ages obtained using the ^^CO (3 — 2) line wings to guide the analysis. Values of 
Pc = 3.7 — 6.2 X km s^^ were obtained for NGC 2264 O and are ~3 times 

smaller than Pqf = 18 x 10^"^ Mq yr^^ km s^^. Overall, the values of Pc and Pqf are similar 
within factors of a few, indicating that the NGC 2264 O outflow has enough force to drive 
the bulk motions of the cloud material. 

Performing a similar analysis for NGC 2264 H results in Pc = 3.7 — 7.3 x 1O~^M0 yr~^ 
km s~^. Compared with the force of the outflow, Pqf = 18 x lO~^M0yr~^ km s"^, Pc is 
several times smaller. Calculations made using line wings of the ^^CO (3 — 2) emission for 
the outflow ages and energetics, find values of Pc — 2.1 — 2.2 x 10~^ Mq yr~^ km s~^, similar 
to the force of the outfiow, Pqf — 2.4 x lO~^M0yr~^ km s"^. Values of Pc and Pqf are 
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within an order of magnitude, similar to NGC 2264 O. Again it appears that the outflow- 
force is capable of driving bulk motions in the cloud. 

4.2. Energetics Comparison 

4-2.1. Turbulent Energy 

The turbulent energy contributed to the cloud by the NGC 2264 O & H outflows, Et^rb, 
was estimated using our ^^CO (3 — 2) data (observations of an optically thin species, such as 
C-'^^O, were not available). The estimate was made by comparing average FWHM velocity 
values at outflow wing emission peaks. Aw (on OF), to ambient cloud regions away from the 
outflow, Av(Ambient). In all regions, the line wings were excluded from the Gaussian flt. 
Since the ^^CO (3 — 2) is optically thick, opacity broadening due to saturation effects should 
be considered (Phillips et al. 1979). Compared to an optically thin species, saturation 
is expected to broaden the FWHM of the ^^CO lines by factors of 2.6 and 2.2 toward 
IRAS 25 and 27, respectively. Normalizing our measured FWHM values by these factors 
and assuming that the ^^CO opacity is to first order constant over the mapped regions, we 
derive values of Af (onOF) = 0.81 kms~^ and Af (Ambient) = 0.54 kms~^ for ^^CO (3 — 2) 
emission surrounding IRAS 25. Normalized FWHM values of At'(onOF) = 0.86 kms~^ and 
Af (Ambient) = 0.45 kms~^ were obtained for IRAS 27. 

FWHM velocity values measured toward peak outflow emission probe many cloud char- 
acteristics including turbulent velocity structures {vfurb), thermal broadening (kT/m), and 
potentially components due to cloud magnetic fleld (vb): Af(onOF) = {v^^^^^ + kT/m + 
v^y^"^- The turbulence associated with the FWHM velocity of the outflowing gas, Vturb, has 
two components. These are general cloud turbulence, Vturb cloud, present both toward and 
away from the outflow, and turbulence due to the outflow itself, Vturb of- These components 
combine quadratically, ■u^^^^ = ■u^^^^ ^^^^^ + ■u^^^^ Qp. The FWHM velocity measured toward 
ambient cloud regions away from outflow emission is Av (Ambient) = {v^^^^ ^^^^^ + '^'s + 
kT/my/'^. Combining these equations and correcting for opacity broadening effects, the 
turbulent velocity due to outflow motions, Vt^rb of-, was obtained from FWHM velocity mea- 
surements: Qp = [Af (onOF)]^ — [Af (Ambient)]^. The turbulent energy contributed to 
the cloud by the molecular outflow was found using Eturb = 1/2 Mqf Vturb of where Mqf 
is the total outflow mass, obtained by summing contributions of red and blue outflow lobes 
(see Table 8). For regions surrounding IRAS 25 and 27, we found Eturb — O.75M0km^s~^ 
and 0.37 Mq km^ s~^, respectively. 

The coupling efficiency of outflow energy to the ambient cloud was estimated by com- 
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paring the turbulent cloud energy contributed by outflows, Eturb, and outflow kinetic energy, 
Eqf (see Table 9). Values of Eqf are orders of magnitude larger than Eturb, resulting in 
coupling efficiencies of < 0.5% in the regions surrounding IRAS 25 and 27. The low coupling 
efficiency of molecular outflow energy to the surrounding environment together with the re- 
sults of our centroid velocity analysis suggests that most of the outflows' energy is directed 
along the flow axis and deposited in the intercloud medium leaving adjacent regions of the 
parent cloud largely undisturbed. 



4-2.2. Gravitational Binding Energy 

Clump gravitational binding energy, Eqc, was obtained using results of the continuum 
analysis (Table 5) and the expression for gravitational potential energy of a uniform gravitat- 
ing sphere of constant density, Eqc = 3/5 (G M^^/R). Here G is the gravitational constant, 
Mff^ is the molecular mass, and R is the sphere radius used to compute Mh^. For IRAS 25, 
the average clump radius is ~ 34" and substitution yields Eqc — 640MQkm^s~^. 

The continuum image of the region surrounding IRAS 27 contains three sources, 27 SI, 
27 S2, and 27 S3. Core masses, Mh2, derived from 870 /xm flux measurements were used to 
individually determine contributions to Eqc for the three sources. Since the core masses, 
Mh2: were derived from fluxes observed within 45" apertures, R — 22.5" was used to obtain 
Eac of individual cores. The results were summed to estimate the total clump gravitational 
potential energy: Eqc = 5.5M0km^s~^. The result is included in Table 9 and is a lower 
limit to PEqc since individual sources lie along a ridge of extended low-intensity 870 /xm 
emission. The mass associated with the low-level emission was not found since the FIR SED 
analysis did not apply. 

Table 9 compares outflow kinetic energy, turbulent energy contributed to the cloud by 
outflows, and cloud clump potential energy for the regions surrounding IRAS 25 and 27. The 
molecular cloud appears to remain gravitationally bound, since turbulent energy injected into 
the cloud by outflows, Eturb, is much smaller than the cloud clump potential energy, Egc- 
The low outflow energy coupling efficiency may explain why, even in the presence of multiple 
outflows, the cloud cores appear to maintain their overall integrity, with Eqc ^ Eturb- The 
dismanthng effects the NGC 2264 O & H outflows are having on the molecular cloud appear 
to be largely localized to the portions of the clouds contained within opening angles of the 
outflow lobes. 
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5. Summciry 

We have made 870 fim continuum measurements and ^^CO (3 — 2), ^^CO (3 — 2), and 
^^CO (1 - 0) observations of IRAS 25 & 27 in the northern cloud complex of NGC 2264 
using the HHT and 12m telescopes. Continuum 870 /xm OTF maps (5'x5') were made 
toward the sources IRAS 25 (IRAS 06382+1017) & 27 (IRAS 06381+1039) and several 
continuum cores were identified. SEDs were constructed for the cores and used to derive 
column densities, gas masses, FIR luminosities, and dust temperatures. A molecular line 
and energetics analysis of the outflows associated with NGC 2264 O & H was performed and 
multiple outflows were found within each region. In order to assess the impact star formation 
is having on the NGC 2264 cloud, the turbulent energy, outflow dynamical energy, and cloud 
clump gravitational potential energy were compared. The results are consistent with the 
interpretation that the cloud complexes are maintaining their overall integrity except along 
outflow axes where cloud material directly interacts with the outflows. The outflows deposit 
most of their energy outside of the molecular cloud, resulting in a weak < 0.5% couphng 
between outflow kinetic energy and cloud turbulent energy. 

We are grateful to C. Kulesa, G. Wolf- Chase, B. Vila-Vilaro, and the Arizona Radio 
Observatory staff and telescope operators for their helpful discussions and aid in assembling 
this data set. This work was supported by NASA GSRP grant number 6267. 
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Table 1. Spectral Line Observations 
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(") 
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(K) 
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Table 2. 870 /im Source Flux 



Source 


ai9B0 


5l950 


45" Aperture Flux 


3' X 5' Aperture Flux 








(Jy) 


(Jy) 


IRAS 25 


06:38:17.0 


10:18:00.0 


2.60±0.52 


8.29±1.66 


25 SI 


06:38:17.8 


10:18:02.6 


6.36±1.27 




25 S2 


06:38:19.2 


10:17:47.7 


5.99±1.20 




IRAS 27 (27 SI) 


06:38:13.0 


10:39:45.0 


3.81±0.64 


22.38±4.48 


27 S2 


06:38:13.0 


10:39:01.0 


2.24±0.45 




27 S3 


06:38:15.9 


10:38:22.5 


1.96±0.39 





Table 3. Source Flux 



Source 


Wavelength 


Flux 


a RMS 


Aperture 


References 




(^im) 


(Jy) 


(Jy) 


(") 




IRAS 25 


870 


8.29 


1.66 


180x300 


1 




850 


6.25 


1.25 


180x300 


2 




450 


21.9 


4.38 


180x300 


2 




100 


86.66 


30.33 


180x300 


3 




60 


43.45 


8.69 


90x285 


3 




25 


7.908 


0.949 


45x276 


3 




12 


1.283 


0.192 


45x270 


3 


IRAS 27 (27 SI) 


1300 


0.8 


0.16 


27.3 


2 




870 


3.81 


0.76 


45.0 


1 




850 


3.20 


0.64 


45.0 


2 




450 


31.8 


6.36 


45.0 


2 




100 


144.9 


37.67 


180x300 


3 




60 


53.66 


11.81 


90x285 


3 




25 


2.755 


0.193 


45x276 


3 




12 


0.25 


0.001 


45x270 


3 


27 S2 


1300 


0.7 


0.14 


27.3 


2 




870 


2.24 


0.45 


45.0 


1 




850 


1.87 


0.37 


45.0 


2 




450 


12.1 


2.42 


45.0 


2 




100 


144.9 


37.67 


180x300 


3 




60 


53.66 


11.81 


90x285 


3 




25 


2.755 


0.193 


45x276 


3 




12 


0.25 


0.001 


45x270 


3 


27 S3 


870 


1.96 


0.39 


45.0 


1 




850 


1.72 


0.34 


45.0 


2 




450 


12.9 


2.58 


45.0 


2 




100 


144.9 


37.67 


180x300 


3 




60 


53.66 


11.81 


90x285 


3 




25 


2.755 


0.193 


45x276 


3 




12 


0.25 


0.001 


45x270 


3 



Note. — 



References: (1) This work. (2) Wolf-Chase et al. (2003) (3) IRAS PSC 
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Table 4. Parameters of FIR SED Model Fits 



Source 


T 

(K) 


''"870 










id 




(Lo) 


IRAS 25 


10 


1.3 X 10-2 


1.4 


1.2 X 10-8 


2.4 X 10" 


-9 


21 X 10- 


-2 


88 




40 


1.5 X 10-3 


1.4 


1.2 X 10-8 


2.0 X 10" 


■10 


1.7 X 10" 


-2 


88 


IRAS 27 (27 SI) 


30 


2.7 X 10^3 


1.4 


1.7 X 10-8 


9.3 X 10- 


-10 


5.3 X 10" 


-2 


120 


27 S2 


36 


1.3 X 10-3 


1.2 


1.7 X 10-8 


3.2 X 10" 


-10 


1.9 X 10" 


-2 


92 


27 S3 


33 


1.2 X 10-3 


1.4 


1.7 X 10-8 


5.1 X 10" 


-10 


2.9 X 10" 


-2 


100 



Table 5. Core Properties 







Source 








Ay 










(Mq) 


(cm-2) 


(cm-3) 










IRAS 25 


140 


3.0 X 10^3 


3.0 X 10^ 


330 








IRAS 27 (27 SI) 


12 


7.0 X 10^2 


8.8 X 10* 


77 








27 S2 


3.4 


1.9 X 10^2 


2.5 X 10^ 


21 








27 S3 


5.5 


3.2 X 1022 


4.0 X 10* 


35 








Table 6. 


C0(3 


— 2) Line Wings 






Source 


F BW 


HV BW 


LV BW 


Line Core 


LV RW 


HV RW 


F RW 




(km s-i) 


(km s-i) 


(km s-i) 


(km s-i) 


(km s-i) 


(km s-i) 


(km s-i) 


IRAS 25 


-10.5, 5.6 


-10.5, -2.5 


-2.5, 5.6 


5.6, 8.8 


8.8, 19.0 


19.0, 27.5 


8.8, 27.5 


IRAS 27 


-20.0, 6.7 


-20.0, 0.0 


0.0, 6.7 


6.7, 9.8 


9.8, 20.0 


20.0, 28.5 


9.8, 28.5 



Table 7. LTE Outflow Analysis Properties 



Source 


Line Wing 


Tex 

(K) 


~32 


-^i-*, thin 

(xlOlS cm-2) 


Nh 

(xl02i cm-2) 


fthin 


IRAS 25 


F RW 


20 


0.16 


14 


8.6 


0.17 


IRAS 25 


F BW 


20 


0.28 


22 


13 


0.10 


IRAS 27 


F RW 


30 


0.10 


14 


8.5 


0.04 


IRAS 27 


F BW 


30 


0.03 


6.6 


3.9 


0.06 



Table 8. NGC 2264 ^^C0{3 - 2) Outflow Energetics 



Outflow 


(kms-l) 


MoF 

(M©) 


POF 
(Mq kms^i) 


EoF 

(xlO^ M© 


LoF 

(L©) 


MoF 
(xlQ-^'Moyr-i) 


^dyn 

(x 10^ yr) 


Pop 

(xlO-*M0yr-l km s'^) 


O (LV RW) 


11 


1.5 


16 


0.9 


0.81 


8.1 


1.8 


8.1 


O (LV BW) 


10.5 


2.9 


31 


1.6 


0.89 


9.8 


3.0 


10 


O (F RW) 


19.5 


1.5 


29 


2.8 


4.5 


14 


1.0 


28 


O (F BW) 


18.5 


2.9 


54 


5.0 


4.9 


17 


1.7 


32 


H (LV RW) 


12 


0.59 


7.1 


0.43 


0.13 


1.1 


5.4 


1.3 


H (LV BW) 


8.0 


0.79 


6.3 


0.25 


0.07 


1.4 


5.7 


1.1 


H (F RW) 


20.5 


0.59 


12 


1.2 


0.65 


1.9 


3.2 


3.8 


H (F BW) 


28 


0.79 


22 


3.1 


3.1 


4.9 


1.6 


14 



Note. — 



Nomenclature: O = NGC 2264 O and H = NGC 2264 H 



-22- 



Table 9. Cloud Core Energetics 



Source 




MoF 


Eac 






Comparison 




(M©) 


(M©) 


(M© Ws-^) 


(M©km2s-2) 


(M©km2s-2) 




IRAS 25 


140 


4.1 


640 


430 - 780 


0.75 


Egc ~ 850 


IRAS 27 


21 


1.4 


5.5 


68 - 250 


0.37 





Note. — BoF values were determined using ^^CO (3 — 2) and ^^CO (3 — 2) line wings, respectively. 
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SED for IRAS 25 
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Fig. 1.— FIR SEDs for sources IRAS 25, IRAS 27 (27 SI), 27 S3, & 27 S2 (clockwise from 
upper left). Prom the shortest wavelength, fluxes included are IRAS 60 pLm and 100 //m data, 
450 //m and 850 //m SCUBA measurements, 870 //m HHT data from this study, and where 
available (IRAS 27 and 27 S2), 1300 jim fluxes. For IRAS 25, a two component BB function 
was flt to this SED. All other sources are well described by single temperature functions. 
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IRAS 25: ^^CO &c '■'C0(3-2) 
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IRAS 27: ^^CO & '^C0(3-2) 
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Fig. 2.— & 13Q0 (3 - 2) overlayed spectra of IRAS 25 and 27. Vertical lines show LV 
and HV red and blue wings and line core (see Table 6). 
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Fig. 3. — Contours of red (solid) and blue (dashed) ''"^00(3 — 2) outflow emission associated 
with NGC 2264 O & H overlayed on 870 /xm continuum maps (thick solid lines). Locations 
of sources within each region are shown, including IRAS 25 & 27 (squares), submillimeter 
continuum objects ("X"), HH 124 emission knots (asterisks), and VLA sources (triangles). 
Continuum 870 /xm contours begin at 2a = 0.38 Jy and 3a = O.lGJy for IRAS 25 & 27, 
respectively, increasing in steps of 0.29 K & 0.27 K toward emission peaks, 1.85 Jy and 2.25 Jy. 
F RW contours begin at 6.6Kkms~^ for NGC 2264 O & H increasing in 2.6Kkms~^ steps 
toward emission peaks, 36Kkms~^ and 31Kkms~^, respectively. F BW contours start at 
6.3 Kkms~^ and 6.8 Kkms~^ for NGC 2264 O & H, respectively, and increase by 1.3 Kkms~^ 
and 2.6Kkms-^ Peak F BW emission is 17Kkms-^ and 32Kkms-i for NGC 2264 O & 
H, respectively. 
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Fig. 4. — Contours of the HV red (solid) & blue (dashed) ^^CO (3 — 2) emission (see Table 6) 
associated with NGC 2264 O & H. HV blue contours begin at 2o" = 1.6 K km s^^ and 
2(7 = 1.5 K km s^^, increasing in steps of 0.15 K km s^^ and 0.7 K km s"^ for NGC 2264 O 
& H, respectively. Contours of HV red emission begin at 2 a = 1.8 K km s~^ and 2 a = 1.0 K 
km s~^, increasing by 0.8 K km and 1.0 K km s^^. HV emission peaks at 2.2 K km s^^ 
and 5.4 K km (blue) and 8.2 K km s-^ and 7.0 K km s'^ (red) for NGC 2264 O & H, 
respectively. 
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Fig. 5.— Model of the ^^C0{3 - 2) outflows surrounding IRAS 25 & 27 (left to right). 
The red (solid) & blue (dashed) lobes of the outflows and their approximate positions with 
respect to half power 870 /im continuum contours are indicated. Two molecular outflows 
are present in each region. The sizes, orientations, and locations of IR reflection nebulae 
associated with IRAS 25 & 27 are indicated with sohd arcs. 
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Fig. 6. — Velocity centroid maps of ^^CO (3 — 2) line core emission associated with 
NGC 2264 O & H. Velocity centroid contours (thin solid lines) are shown and labeled (km 
s~^). Half power contours of red and blue outflow emission are illustrated with thick solid 
and dashed contours, respectively. 



